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Abstract 

We update the theory predictions for the mass difference AM S , the width difference Ar s and the 
CP asymmetry in flavour- specific decays, a| s , for the B s — B s system. In particular we present 
a new expression for the element T s 12 of the decay matrix, which enters the predictions of Ar s 
and af s . To this end we introduce a new operator basis, which reduces the troublesome sizes of 
the 1 jm h and a s corrections and diminishes the hadronic uncertainty in Ar s /AM S considerably. 
Logarithms of the charm quark mass are summed to all orders. We find Ar s /AM S = (49.7 ± 
9.4) ■ 10" 4 and Ar s = (/ Ba /240 MeV) 2 [(0.105 ± 0.016) 5 + (0.024 ± 0.004) B' s - 0.027 ± 
0.015] ps _1 in terms of the bag parameters B, B' s in the NDR scheme and the decay constant fs s ■ 
The improved result for Y\ 2 also permits the extraction of the CP- violating B s —B s mixing phase 
from af s with better accuracy. We show how the measurements of AM,, AT S , af s , Aqp(B s — > 
J/ijj(p) and other observables can be efficiently combined to constrain new physics. Applying 
our new formulae to data from the D0 experiment, we find a 2a deviation of the B s —B s mixing 
phase from its Standard Model value. We also briefly update the theory predictions for the 
B d -B d system and find Ar d /AM d = (52.6^2^) ■ 10~ 4 and af s = (-4.8±};°) ■ 10~ 4 in the 
Standard Model. 



PACS numbers: 12.38.Bx, 13.25.Hw, 11.30Er, 12.60.-i 
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1 Introduction 

Flavour-changing neutral current (FCNC) processes are highly sensitive to new physics around 
the TeV scale. Global fits to the unitarity triangle show an excellent agreement of b — > d and s — > 
d transitions with the predictions of the Cabibbo-Kobayashi-Maskawa (CKM) mechanism [1,2]. 
Extensions of the Standard Model can contain sources of flavour-changing transitions beyond the 
CKM matrix. Models without these new sources are termed to respect minimal flavour violation 
(MFV). Despite of the success of the MFV hypothesis in b — > d and s — > d transitions there 
is still sizable room for non-MFV contribution in b — > s transitions. For instance, an extra 
contribution to b — > sqq, q — u,d, s, decay amplitudes with a CP phase different from &rg(V t *Vtb) 
can alleviate the ~ 2.6a discrepancy between the measured mixing-induced CP asymmetries in 
these b — > s penguin modes and the Standard Model prediction [3]. Models of supersymmetric 
grand unification can naturally accommodate new contributions to b — ► s transitions [4]: right- 
handed quarks reside in the same quintuplets of SU(5) as left-handed neutrinos, so that the large 
atmospheric neutrino mixing angle could well affect squark-gluino mediated b — > s transitions 
[5]. 

Clearly, B s — B s mixing plays a preeminent role in the search for new physics in b — ► s 
FCNC's. B s — B s oscillations are governed by a Schrodinger equation 




W)) \ _ f M s_i r s\ ( \B a (t)) 

\B.(t)) i \r 2 1 ) { \B s (t)) 



(1) 



with the mass matrix M s and the decay matrix V s . The physical eigenstates \B H ) and \B L ) with 
the masses M H , M L and the decay rates T H , T L are obtained by diagonalizing M s — iT s /2. The 
B s —B s oscillations in Eq. (Q]) involve the three physical quantities | Aff 2 1 , |rf 2 | and the CP phase 
<p s = arg(— Mf 2 /r* 2 ) (see e.g. [6]). The mass and width differences between Bl and Bh are 
related to them as 

AM, = M S H — M[ = 2\M? 2 \, AT S = T S L — T S H = 2|r s 12 |cos0 s , (2) 

up to numerically irrelevant corrections of order m^/M^. AM S simply equals the frequency 
of the B s — B s oscillations. A third quantity providing independent information on the mixing 
problem in Eq. (Q~|) is 

c -L i q 1 1 i o I ^ s 

at = Im = sin 6 S = tan 6 S . (3) 

fs Mf 2 |Mf 2 | r AM S v v 

af s is the CP asymmetry in flavour-specific B s — >• / decays, which means that the decays B s — > / 
and B s — >• / (with / denoting the CP-conjugate final state) are forbidden [7]. The standard way 
to access a| s uses B s — > X s £ + T7^ decays, which justifies the name semileptonic CP asymmetry 
for af s . (See e.g. [6, 8] for more details on the phenomenology of B s — B s mixing.) 

It is important to note that new physics can significantly affect Mf 2 , but not rj 2 , which is 
dominated by the CKM-favoured b — > ccs tree-level decays. Hence all possible effects of new 
physics can be parameterised by two real parameters only, for instance |Mf 2 | and 4> s . While 
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b u,c,t s 




s u,c,t b 



Figure 1: In the lowest order Mf 2 is calculated from the dispersive parts of the box diagrams 
on the left. It is dominated by the top contribution. The result involves only one local \AB\ = 
2 operator, shown in the right picture. The leading contribution to T s 12 is obtained from the 
absorptive parts of the box diagrams on the left, to which only diagrams without top quark line 
contribute. To lowest order in the heavy quark expansion two \AB\ = 2 operators occur, the 
A/rrib corrections involve five more. 



|Mf 2 | is directly related to AM S , the extraction of <j) s from either Ar s or a| s requires an accurate 
knowledge of rf 2 . 

In the Standard Model M( 2 and TJ 2 are computed from the box diagrams in Fig. [Hand QCD 
corrections in the desired order. The Standard Model prediction for M i2 reads: 

(~<2 tut 

M 12 = ^f^M^(V tb V t :) 2 f} B S {x t )f B B, (4) 

where Gf is the Fermi constant, the V*/s are CKM elements, M Ba and Mw are the masses of 
B s meson and W boson and the short-distance information is contained in fj B So(x t ): So(x t ) is 
the Inami-Lim function, which depends on the top mass m t through x t = mf/M^, and r/s is a 
numerical factor containing the leading and next-to-leading QCD corrections [9]. The calculation 
of M 12 involves the four-quark operator (a, (5 = 1, 2, 3 are colour indices): 

Q = ^7^(1 - 7 5 )6 Q s /3 7 / "(l - 7 5 )&/3. (5) 

All long-distance QCD effects are contained in the hadronic matrix element of Q and are param- 
eterised by f Bs B'. 

(B S \Q\B S ) = *-Mlj 2 B B. (6) 

The recent observation of the B s — B s mixing frequency AM S = 2|Mf 2 | at the Tevatron [10] 
yields a powerful constraint on extensions of the Standard Model [11-14]. The results from the 
D0 and CDF experiments obtained with 1 fb 1 of data, are [15] 

17ps _1 < AM a < 21PS" 1 @90%CL D0 

AM S = 17.77 ±0.10 (syst) ± 0.07 (gM) ps _1 CDF. (7) 
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Figure 2: Leading-order CKM-favoured contribution to T s l2 , which arises from B s decays to final 
states (indicated by the dashed lines) with a (c, c) pair and zero strangeness. The crosses denote 
any of the operators Q\-% of the | AB\ = 1 hamiltonian. The Cabibbo-suppressed contributions 
correspond to diagrams with one or both c quarks replaced by u quarks. 



While the precise measurement in Eq. © sharply determines |Mf 2 |, the uncertainty of /§ B, 
which is around 30%, blurs the extraction of some new physics contribution adding to So(x t ) in 
Eq. (HJ. Alternatively one can study the ratio AMd/AM s , where AMj is the mass difference 
in the Bd — B d system. While the hadronic uncertainty in the ratio f^Bj (fg d BB d ) is smaller, 
one is now dependent on | V^/V^l 2 . Even if one assumes non-standard contributions only in B s 
physics, but not in the quantities entering the global fit of the unitarity triangle, | Vjd/V^| 2 is only 
known to roughly 40% [2] leaving equally much room for new physics in |Mf 2 |. 

Adding experimental information from AY S or af s helps in two ways; first, one can study 
the CP-violating phase <p s , which is totally unconstrained by AM S , through Eqs. © and ([3]). 
Second, one expects cancellations of hadronic parameters in the ratio r^ 2 /Mf 2 , which enters 
af s and AT S / AM S . All decays into final states with zero strangeness contribute to rj 2 , which 
is dominated by the CKM-favoured b — > ccs tree-level contribution. In the first step of the 
calculation the W-boson is integrated out and the W- mediated | AB\ = 1 transitions are described 
by the usual effective \AB\ = 1 hamiltonian with the current-current operators Qi, Q2 and the 
penguin operators Q3-Q, Qs [16]. The leading contribution to T\ 2 in this effective \AB\ = 1 
theory is shown in Fig. [21 In the second step one uses an operator product expansion (OPE), the 
Heavy Quark Expansion (HQE), to express I^ 2 as an expansion in the two parameters A/m^ and 
a s (mb). Here a s is the QCD coupling constant and A is the appropriate hadronic scale, which 
quantifies the size of the hadronic matrix elements. The HQE links the diagrams of Fig. [2] to the 
matrix elements of local AB = 2 operators. In addition to the operator Q in Eq. © one also 
encounters 

Qs = s a (l + 75)^5,3(1 + 75)6/3, (8) 

whose matrix element is parameterised by a bag parameter Bs in analogy to Eq. ©. The leading 
contribution to T^ 2 was obtained in [7, 17]. Today rj 2 is known to next-to-leading-order (NLO) 
in both A/m b [18] and a s {m b ) [19,20]. The 1998 result [19] 

(tt) = (2Tntev) 2 ' - 006B+0 - 150 ^-°- 063 i (9 > 
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with the average total width T s = (T S L + T S H ) /2 is pathological in several respects: first, the A/m& 
correction -0.063 is unnaturally large and amounts to around 40% of the total result. Second, 
the coefficient of B cancels almost completely, the result is therefore dominated by the term 
proportional to B s ~ 0.9, so that the cancellation of hadronic quantities from the ratio AT S /AM S 
is very imperfect. Third, both the A/m b and a s corrections, which diminish the coefficient of 
Bs from 0.22 to 0.15, are negative, and these numerical cancellations between leading-order 
(LO) order result and corrections increase the relative uncertainty of the prediction for Ar s /T s . 
In the following section we argue that these pathologies are caused by a poor choice of the 
operator basis used in [18-20] and propose a different basis. We also improve the prediction of 
Ar s /AM S and Ar s /T s in several other aspects, by summing logarithms of the charm mass to 
all orders in a s , by using different renormalisation schemes for the 6-quark mass, by including 
CKM-suppressed contributions and by modifying the normalisation related to the factor 1 /T s 
in Eq. ©. In Sect. [3] we present numerical updates first of AM S , Ar s and a| and then of 
the corresponding quantities in the ^-system. In Sect. @] we show how the expressions for 
the mixing quantities change in the presence of new physics. Here we discuss how to combine 
different present and future measurements to constrain | Mf 2 1 and (f> s and advocate a novel method 
to display the constraints on possible new short-distance physics in B s —B s mixing. Sect. Ogives 
a road map for future measurements and calculations and Sect. [6] summarises our results. 



2 Improved prediction of T 
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We write r? 2 as [21] 



1 12 



\2 -pec I o \ \ I "uc i \2 i -km 
A - 1 12 + Z A c A u J- 12 + A u 1 12 



\ 2 pec 
A t 1 12 



12 



2 A t a u (r 



12 



1 12) 



+ A l (r 



cc 
12 



12 , 



(10) 

(11) 



with the CKM factors Aj = V^Va for i — u,c, t. In Eq. (fTTI) we have eliminated A c in favour of 



A t using X u + A c + A t = to prepare for the study of Y\ 2 /Ml 2 . Since \X U \ <s^ \X t \ 
clearly dominates T^ 2 . For ab = cc, uc, uu we write [19, 21] 



A c ,r 



cc 
12 



-pab 
1 12 



24ttM Bs 



G ab (B s \Q\B s ) - Gf(B s \Q s \B s ) 



+ r 



12,1/mj, 



(12) 



The coefficients G and G a s are further decomposed as 



G 



ab 



pab _|_ p 



ab 



G1 



_pab pab 

S S 



(13) 



Here F and F§ are the contributions from the current-current operators Q 12 while the small 
coefficients P ab and Pg b stem from the penguin operators Q3-6 and Q$. (Note that in [19], where 
only the dominant Tf 2 was considered, these coefficients had no superscript 'cc'.) Numerical 
cancellations render F cc small with |F CC /F| C | ^0.03 which explains the small coefficient of B 
in Eq. ©. 
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We parameterise the matrix element of Qs as 

(B S \Q S \B S ) = -hl 2 B J 2 B B' s . (14) 

Formulae for physical quantities are more compact when expressed in terms of B' s rather than 
the conventionally used bag parameter B s . The two parameters are related as 



M 2 

B's = ,_ *, 2 B S . (15) 



In the vacuum insertion approximation (VIA) the bag factors B and B$ are equal to one. Through- 
out this paper we use the MS scheme as defined in [19,21] for all operators. Therefore the masses 
fn b and m s appearing in Eq. (fT5l) correspond to the MS scheme as well. 

^12 i/ mb comprises effects suppressed by A/m b . We will discuss it later, after transforming 
to our new operator basis. 

2.1 New operator basis 

When calculating T 12 to leading order in A/m b , one first encounters a third operator Qs in 
addition to Q and Qs defined in Eqs. © and ([8]): 

Qs = s a (l + 75)6,35,3(1 + 7 5 )& a , (16) 

However, a certain linear combination of Q, Qs and Qs is a l/m;,-suppressed operator [18]. 
This l/m b -suppressed operator reads 

R = Q s + ocxQs + ^a 2 Q, (17) 
where a 12 contain NLO corrections, which are specific to the MS scheme used by us [19]: 

ax = 1 + — : C/ 12m 1-6 , a 2 = 1 + — ^ Ct 6 In h — . (18) 

47T \ mb / An \ m b 2 / 

Here C/ = 4/3 is a colour factor and /i 2 is the scale at which the operators in Eq. (fT7l) are defined. 
The coefficients G and Gs in Eq. (fT2)) depend on fi 2 and this dependence cancels with the /x 2 - 
dependence of (-B S |Q(// 2 )|-B S ) and | Qs (a* 2) | -S*) - I n lattice computations the /^-dependence 
enters in the lattice-continuum matching of these matrix elements. In our numerics we will 
always quote the results for /i 2 = m b . In [18-20] Eq. (TT7T) has been used to eliminate Qs in 
favour of Ro leading to the result in Eq. ©. The matrix element of Qs reads 

(B s \Qs(^)\B s ) = \mIJIB' s {im 2 ). (19) 



In analogy to Eq. (TT31) we define 

^ = ■_ _ 

(m b (/i 2 ) + m s (n 2 )) 



Ml 



B 'sM = 7=r7^ ,',,~ 2 BsM. (20) 



6 



Theoretical update of B s — B s mixing 



For clarity we have explicitly shown the /^-dependence in Eqs. (1191) and (1201) . which was skipped 
in Eqs. ©,CE]> and (Q3]). In VIA B s = 1 and (B S \Q S \B S ) is much smaller than (B S \Q\B S ) 
and (B s \Qs\B s ). The small coefficient 1/3 in Eq. (fl~9l) is the consequence of a cancellation 
between the leading term in the 1/N C expansion, where N c = 3 is the number of colours, and 
the factorisable 1/N C corrections: 1/3 = 1 — 2/N c . One naturally expects that the bag factor 
Bs substantially deviates from 1. However, a lattice computation found Bs = 0.91 ± 0.08 [22], 
showing that the matrix element of Qs is indeed small. Thus (B S \R \B S ) = kjm b implies a 
strong numerical relationship between B and B s which can be used to constrain B s /B entering 
Ar s /AM S . Yet it is more straightforward to use Eq. (flTT) to eliminate Qs altogether from Ti 2 in 
favour of Qs- The coefficient of B will change and and the coefficient of B' s is expected to be 
small in view of the factor of 1/3 in Eq. (fT9l ). Using further the bag parameters of Eqs. © and 
COD), Yf 2 of Eq. CL2]) now reads 



-pab 
1 12 



G 2 m 2 
24vr BsJBs 



Qob _|_ ^2 (~yab 



\B + Gfa, 



+ T* 1/m . (21) 



The new 1 /m^-corrections are related to \/ mb appearing in Eq. (fl~2l) as 



-pab 

1 12,1/mt 



-pab i G\m\ 

12,l/m b f 247rM Bs 



n ab,(0) 



(S s |i2 |-B s ). 



(22) 



Here we have taken into account that the result of [19, 20] includes the A/m^ terms without 
penguin contributions and to LO in a s : consequently we have changed —G a s h to F^ 6 '^ '', which is 
the LO approximation to Fg b . Recalling \G ab \ \G s b \ and B, B' s rs 1 one easily verifies from 

in 



ab 

12,l/m b 



Eq. (12TI) that the first term proportional to B dominates over the second term. Since T 
Eq. © is negative and the shift in Eq. (|22|) adds a positive term our change of basis also leads 
to |r^ jy m J < |r^ x/ m J- Further the a s -corrections contained in ai 2 , which multiply G^ 6 '^ ^ 
in Eq. (|2TI) . temper the large NLO corrections of the old result. These three effects combine 
to reduce the hadronic uncertainty in Ar s /AM S substantially. In other words: the uncertainty 
quoted in [19, 20] is not intrinsic to Ar s /AM S but an artifact of a poorly chosen operator basis. 



2.2 A closer look at 1/ nib corrections 

At order l/m b one encounters the operators R of Eq. (flTT) . 



Ri = 


m s 






m b 






1 




R 2 = 




s 




m 2 b 






1 




R 3 = 




s 




m\ 





(23) 



and the operators Ri which are obtained from the i?j's by interchanging the colour indices a and 
(3 of the two s fields [18]. At order l/m& only five of these operators are independent because of 
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relations like R 2 = —R2 + 0(1/ ml). Writing (for ab = cc, uc, uu) 

>(B s \Ro\B s ) + j^[gf{B s \R j \B s )+gf{B s \R j \B s ) 



12,l/m 6 



G 2 F ml 
24ttM Bs 



9o bl 



(24) 



the coefficients gf and gf read [18,21,23]: 



gl c = VT^Tz(l + 2z)C ( 2 m + F c s c{0) = v / r 3 4i(l + 2^)C| 0) 3Cj 0) + 2G 



?cc(p) 



(0) 



.(0) 



9i 



-2 V / T^4i(l + 2z) C[ 0) 3Cf ) + 2C\ 



,(0) 



.(o) 



-2 



cr 3cr + 2c 2 u 



9t 



\-2z- 2z 2 r ,,o ) o/ ., (U) o/ ., (0; 
-24 



(o) 



9T 



l-z) 2 (l + 2z)Cf 2 



_|_ jpuc(0) 



9T 



9? 



9 s 



[l-z 



-2V1 - 4*(1 + 2z)C. 
2 1 - 2^-2^ 2 ^(o) 2 



(0)2 



v/r^4i 

v2 



-24 



2T 



^(0)2 



(25) 



. 2 (l + 2^)Cf ) [3Cf + 2Cf 



< c = -2(1 - 


-z) 2 (l + 2^)Cf 


~3C{ 0) 


+ 2Cf" 


9T = 


-2(1- 


<?r = -2(1- 


(l + z + z 2 )C{ 0) 


3Cf 


+ 2^ or 


9? = 


-2(1- 




-12 (l-z) z 2 C{ 0) 


3Cf 


+ 2C? y 


9T = 


-12(1 



z) 2 a 



2z)Cf 2 



,(0)2 



z 2 )Cf 2 
(26) 



and g 



uu 
3 



gf(z 



0) 



gf(z 



0). Here 



m| _ [m e (m e )] 2 
m 2 [m b (m b )] 2 



(27) 



and C^ - 1 ~ —0.3 and C 2 ' ~ 1.1 are the LO Wilson coefficients of the AB = 1 operators Qi 
and a. |16|. 

The contributions involving Ri, Ri, R 3 and R 3 are suppressed by powers of m s /m b or z 2 
and are numerically negligible. The only two important l/m b operators are R and R 2 = —R2 + 
0(1 /ml). As a consequence of the elimination of Q s in favour of Q s no term involving the large 
coefficient C 2 ^ 2 occurs in g$ b . The contribution from R is substantially diminished, and this 
can be understood in terms of a systematic expansion in 1/N C : the coefficients g$ b are colour- 
suppressed due to C\ ~ 1/N C , while they were colour-favoured in the old basis. Since radiative 
corrections cannot change the colour counting, this feature must persist in the yet uncalculated 
order a s /m b . In other words, by changing to our new basis we have absorbed the corrections of 
order N®/m b into the leading order of the l/m b expansion. This improves our result over the 
one in the old basis by a term of order N c a s /m b . (Recall that a s ~ 1/N C , so that N c a s /m b ~ 
N®/m b .) This term (which constitutes a parametrically enhanced correction) would appear, if the 
calculation of a s /m b were done in the old basis. In fact, this term occurs in the NLO calculation 
of [19-21] in the coefficient of Q s but is dropped once Q s is traded for R , because all a s /m b 
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terms are consistently discarded. With the use of our new basis no corrections of order N c a s /rrib 
to g^ b can occur. This feature can also be understood by realising that the large-iV c contribution 
to stems from the right diagram in Fig. [2] with two insertions of Q 2 plus additional planar 
graphs with extra gluons. These diagrams contribute to the coefficients of Q and Qs, but not to 
the coefficient of Qs- (This is easy to see, if one inserts the two Q 2 ' s in the Fierz-rearranged 
form.) Upon elimination of Qs in favour of Rq, the color-suppressed coefficient g ab of Qs 
becomes the coefficient of R . At order \/m b one has to include the momentum of the s quark in 
that diagram and finds a contribution to the g 4 ab 's at order N®. These terms are identical in both 
bases. Our numerical analysis in Sect. [3] follows the pattern revealed by the 1/N C expansion, 
finding the numerical relevance of R drastically reduced compared to the old basis, so that the 
only remaining important l/m& operator is R 2 - 

In the new basis the l/m b corrections have their natural size of order A/m b ~ 20%. To 
be conservative, we have estimated the I /ml terms to verify that this result is not accidental. 
We have found two types of contributions: the first type is calculated by expanding the results 
of Fig. [2] to the next order of the s-quark momentum, yielding operators with more derivatives 
acting on the s quark field. We find that these contributions have the same suppression pattern as 
the gf"s and gf"s. The second type of l/rrc| operators involve the QCD field strength tensor 
and has no counterparts at lower orders. We find small coefficients here as well. Since the size of 
the l/m 2 corrections is well below the uncertainty which we obtain by varying the bag factors of 
the operators in Eq. (|23l ), there is no reason to include these corrections into our numerical code. 

We parameterise the matrix elements (R4) = (B 8 \Ri\Bg)'s as 



(Ro) 



4 
3 



Ml 



m pow2 ^ + m J mb ^ 



- 1 



6 m b 



(R2) = 

(Rs) 
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3 
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Ml 

pow 2 
m b 

Ml 



- 1 
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pow 2 



m bJb s b R^ 
MbJbsBrsi 



m bJb s Bro, 



(Ri) 

(R2) 



- —M 2 f 2 B~ 

3 m b BjBa 



2 
3 

<«,> - I 



M 2 



B 3 



pow 2 
m b 

M 2 s 

pow 2 



- 1 



MlJl B ^ 
MlflB,. (28) 



As usual the bag parameters Br , . . . , B^ parameterise the deviation of the matrix elements 
from their VIA results derived in [18]. The numerical values of the (Ri)'s depend sensitively on 
the choice of the mass parameter ml° w in Eq. (l28l . Clearly, m£ ow is a redundant parameter, as 
any change in m^ ow can be absorbed into the bag parameters. It merely serves to calibrate the 
overall size of the l/m^-suppressed matrix elements such that the bag factors are close to 1. A 
future NLO calculation of the coefficients in Eq. (|26l) will allow us to replace ml° w by a well- 
defined (i.e. properly infrared- subtracted) b pole mass. Our numerical value for ml° w is guided 
by the requirement that the terms in square brackets in Eq. ([28]) are of order 2A/m^ ow ~ 0.2, 
which leads to the estimate ml° w fa 4.8 GeV. A better justification can be given by noting that 
the lattice computations of B, B s and B s in [22] allow for an estimate of (R ) (which may 
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become a determination, once the lattice-continuum matching of (R ) is done at NLO): 

l / r, 1 

B Ro = —B'v + a 2 B--B' Q 1 



f B' s + a 2 B - \ B ' S 



Ml 



ml° w2 (l + m s /m b ) 2 



(29) 



With the central values for B, B,s and B,s given in [22] and the choice m£ ow = 4.8 GeV one finds 
Br = 1.1, while those of the new preliminary lattice computation of [24] imply Br = 1.7. 
Our quoted numerical results in Sect. [3] correspond to conservative ranges for both m^ ow and the 
B R /s. We note that the only places where we use m po> 
is not used in the overall factor mf of \/ mb 
analysis in [21]. 



; , are the matrix elements in Eq. (|28l ; it 
in Eq. (l24l . This is a change compared to the 



2.3 Summing terms of order a n z hi" z 

° s 

The coefficients G ab and Gg b in Eq. ([5]) depend on quark masses through z defined in Eq. (|27l) . At 
order a n s the dominant ^-dependent terms are of the form a n s z ln n z. In [25] and [21] it has been 
shown that these terms are summed to all orders n = 1, 2, . . ., if one switches to a renormalisation 
scheme which uses 

z = (30) 

[m b (m b )} 

Since z is roughly half as big as z, this also reduces the dependence of the coefficients on the 
charm mass. We illustrate the effect for Yf 2 with a numerical example: In the two renormalisation 
schemes one finds 

= (3.3 - 11.4 z + 1.5 z hiz) ■ 10" 3 ps _1 + O (z 2 ) 
Yf 2 = (3.3 - 11.42?) ■ lO^ps" 1 + O (i 2 ) . (31) 

The numerical input is taken from Eqs. (13214381 and Eq. (|39l) below. From Eq. (1311) one verifies 
that the use of z eliminates the z In z term. This issue is particularly relevant for af s and af s , 
which are of order z. The final numbers for all quantities quoted below involve z. We only revert 
to a scheme using z to compare with the previously published results in [19, 20]. 



3 Numerical predictions 

3.1 Input 

For the numerical analysis we use the following set of input parameters: The quark masses 
are [26] 

m b (m b ) = 4.22 ±0.08 GeV m£ ole = 4.63 ± 0.09 GeV (32) 

m pow = 4.8+0.0 GeV 
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tti 2 (in ) 

m c (m c ) = 1.30 ± 0.05 GeV => z = _l _ c = 0.095 ± 0.008 , (33) 



m 2 b {m b ) 
m 2 ijrib) 



0.048 ± 0.004 



ml(m b ) 

m s (2GeV) = 0.10 ±0.02 GeV m s (m b ) = 0.085 ± 0.017 GeV 

m poie = 171,4 ±2.1 GeV => m t (m t ) = 163.8 ± 2.0 GeV (34) 

We will need the meson masses [27] 

M Bd = 5.279 GeV , M Ba = 5.368 GeV . (35) 

The average width T s of the B s mass eigenstates is computed from the well-measured B d life- 
time, 

r Bd = 1-530 ±0.009 ps, (36) 
using T s = l/r Bd (1.00 ± 0.01). Our input of the CKM elements is [2] 

\V US \ = 0.2248 ±0.0016, |V c6 | = (41.5 ± 1.0) • 10" 3 
] "'" = 0.10 ±0.02, 7 = 1.05lg;i2. (37) 



For all predictions within the standard model we assume unitarity of the CKM matrix and we 
determine all CKM elements from the four parameters |14 s |, \V c b\, \V u b/V c b\ and 7. The W 
mass [27] and the strong coupling constant are [28] 

M w = 80.4 GeV , a s (M z ) = 0.1189 ± 0.0010 . (38) 

We note that in the B s system CKM parameters other than \V c b\ (which basically determines 
I V ts |) play a minor role. The same is true for the strange quark mass in Eq. (|34l ). 

The dominant theoretical uncertainties, however, stem from the non-perturbative parameters 
discussed below and from the dependence on the unphysical renormalisation scale We use 
the central values fii = /x 2 = nib and we vary /ii between m^/2 and 2m&. The dependence on /i 2 
is related to the determination of the hadronic quantities and uncertainties associated with fx 2 are 
contained in the quoted ranges for these quantities. 

The situation of the non-perturbative parameters - the decay constant and the bag parameters - 
is not yet settled. Different non-perturbative methods result in quite different numerical results. 
QCD sum rule estimates were obtained for the decay constant f Bs [29], for the bag parameter B 
[30,31] and for Bs [31]. The same quantities have been determined in quenched approximation 
in numerous lattice simulations, see [32] for a review. The only determination of B s was done in 
a quenched lattice simulation in [22]. Unquenched (rif = 2) values are available for f Bs [33,34], 
for B [34, 35] and for Bs [35, 36]. For the decay constant f Bs even a lattice simulation with 2+1 
dynamical fermions is available [37]. 

Unfortunately it turns out that the predictions for f Bs vary over a wide range, 0(200 ± 20 MeV) 
for quenched results, 0(230 ± 20 MeV) for n f = 2, 0(245 ± 20 MeV) for sum rule estimates 
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and (9(260 ± 29 MeV) for rif — 2 + 1, see e.g. [32]. This discrepancy has to be resolved, since 
AM and Ar s depend quadratically on the decay constant! Recently the combinations fs 3 B, 
fs s Bs and f Ba B s were determined for 2+1 light flavors [24]. The authors of [24] claim that the 
combined determination results in a considerable reduction of the theoretical error. 
We will use in our numerics two sets of non-perturbative parameters: 

Set I consists of a conservative estimate for fs s combined with the unquenched determination 
for B [34] and B s [36] and the only published lattice determination of B s [22]: 

f Bs = 240 ±40 MeV 
B = 0.85 ±0.06 => 

B s = 0.86 ±0.08 => 

B s = 0.91 ±0.08 => 

Set II consists of the preliminary determination with 2+1 flavors [24]: 

f Bs ^B = 0.227(17) GeV 
f Ba ^B 1 s = 0.295(22) GeV 

f Ba JW s = 0.305(23) GeV (40) 

The central values of both sets are quite similar, while the errors of set II are smaller by almost a 
factor 3. 

For both sets the bag parameters of the l/m^-corrections are estimated within vacuum insertion 
approximation and we use the following conservative error estimate 

5^ = 1 ±0.5. (41) 

In our computer programs we carefully extract all terms of order o? s and a s /mb, which belong to 
yet uncalculated orders of the perturbation series, and discard them consistently. 



f Bs VB = 0.221(46) GeV 

B' s = 1.34 ±0.12 =^ fB s \jB T s = 0.277(57) GeV 
B' s = 1.41 ±0.12 fB s \l~B~ s ~ = 0.285(60) GeV (39) 



3.2 AM S within the SM 

In the standard model expression (Eq.© & Eq.©) for the mass difference in the £? s -system a 
product of perturbative corrections 07b So) and non-perturbative corrections (f B B) arises. Using 
the above input the perturbative corrections are given by [9] 

= 0.837 (NDR), (42) 

So(xt) = S o r-^\ = AXt ~ ^ ± %l - = 2 ,3 2 7 ± 0.044 (43) 

\ Mfy J 4(1 -x t ) 2 2{l-x t ) i 

Our final values for the standard model prediction 

AM S = (19.30 ± 6.68) ps" 1 (Set I) (44) 
AM S = (20.31 ± 3.25) ps" 1 (Set II) (45) 
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are bigger than the experimental result, but consistent within the errors. Using fs s = 230 MeV 
and the bag parameter from set I, one exactly reproduces the experimental value of AM S . 
The overall error is made up from the following components: 



Input 


AM S 


AM S 




Set I 


Set II 


Ib s 


i +0.361 
1 -0.306 




B 


1 ±0.071 




fl s B 


1 ±0.341 


1± 0.150 


v cb 


1 +U.049 
■•—0.048 


-i +U.049 
J— 0.048 


a s (M z ) 


1 ±0.020 


1± 0.020 


m t 


1 ±0.018 


1 ±0.018 


7 


i +U.0U5 
J— 0.015 


i +U.0U5 
J— 0.015 




1± 0.005 


1± 0.005 


\/E5 2 


1 ± 0.346 


1 ±0.160 



When combining different errors we first symmetrised the individual errors and added them 
quadratically afterwards. The by far dominant contribution to the error comes from the non- 
perturbative parameter f% B. Clearly, in view of the precise measurement in Eq. © it is highly 
desirable to understand the hadronic QCD effects with a much higher precision than today. 



3.3 Ar s , Aiy AM S and af s within the SM 

The main result of this paper is a new, more precise determination of Ti 2 , which is then used to 
determine AT S , AT S / AM S and af s . 

In order to illustrate our progress, we first present the results in the old operator basis used in 
[19,20]. Using the scheme involving m£ ole and z as in [19,20], but updating the input parameters 



to our values in Eqs. 



Ar 



pole 
s,old 



a -ppole,LO 
m s,old 



pole,s 
l fs,old 



( AT, 

\AA'L 



pole 



s / old 



+ 



, we find 



B, 



[0.0025 + 0.0945^- 



240 MeV / 

(0.0335^ + 0.0195 % + 0.0055* 

— ^ — \ [0.0055 + 0.1455' - 
240 MeV/ L s 

(0.0335^ 2 + 0.0195^ + 0.0055 R 



ps" 



ps" 



B' Q 



B 



B 



10.8 + 1.9-f + 0.8-+ 



5 
B' 



B 



0.10 -0.01-f + 0.29 



5 



R 



B 



B' ( 
0.9 + 40.9-^ - (11. I 



5 
5 



Im 



R 2 
B 



+ 8.5 



■ 10" 



5 



2.1- 



Br 
B 



■ 10" 



(46) 
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For simplicity we do not show the uncertainties of the numerical coefficients appearing in the 
square brackets here and in following similar occasions. We assess these uncertainties, however, 
when quoting final results. 

Several comments are in order: in the old basis the coefficient of B in the prediction of 
Ar s is negligible due to a cancellation among AB = 1 Wilson coefficients, thus the term with 
B' s dominates the overall result. This leads to the undesirable fact that the only coefficient in 
AT S / AM S that is free from non-perturbative uncertainties is numerically negligible. Moreover 
in Ar s all 1 /m b -corrections have the same size and add up to an unexpectedly large correction 
(30% of the LO value, 45% of the NLO value). In Eq. (|46l we have singled out the bag factors 
of the two most important sub-dominant operators R 2 and R , while the bag parameters of the 
remaining operators are chosen equal and are denoted by Br. Finally in the old operator basis 
the calculated NLO QCD corrections are large and reduce the final number by about 35% of the 
LO value. 

a% does not suffer from this shortcomings. Here the coefficient without non-perturbative uncer- 
tainties is numerically dominant and the size of the l/m& corrections seems to be reasonable. 
Moreover in this case i? 3 and i? 3 are the dominant subleading operators. Since the overall contri- 
bution of the l/m fo -corrections is relatively small, we choose all bag factors of power suppressed 
operators equal to B R . 

Using the non-perturbative parameters from set I we obtain the following number for Ar s : 

Ar 

Ar s = (0.070 ± 0.042) ps" 1 = AT S ■ r Bd = 0.107 ± 0.065 (47) 

J- s 

This number is in agreement with previous estimates [19,38-40] where different input parameters 
- in particular different values for the decay constant and the bag parameters - were used. In the 
following table we quote the central values of these old predictions and in addition give the 
corresponding results adjusted to the new non-perturbative parameters of set I: 



Reference 


predicted Ar s /r s 


used f Bs 


used B' s 


Ar s /r s (/ Bs = 240MeV, 
B' s = 1.34) 


[19] 


0.054 


210 


1.02 


0.117 


[38] 


0.093 


230 


1.25 


0.114 


[39] 


0.124 


245 


1.36 


0.116 


[40] 


0.118 


245 


1.31 


0.117 



The values in the last column are still bigger than the new number in Eq.(l47]) by about 8%. 
Besides some differences from other input parameters — like quark masses and CKM parame- 
ters — this small overestimate in the last column originates from the use of different methods 
to determine T s in the ratio Ar s /r s compared to this work. Since now very precise values of 
the b-lifetimes are available, we directly use them as an input to determine the total decay rate: 
T s = l/re d . In [19,38-40] we expressed the total decay rate in terms of the semileptonic de- 
cay rate: T s = I^} 16017 / B^ P ■ Doing so (with the 1998 value of B^ p ) one obtains values for 
t b ~ 1.66 ps, which are about 8% larger than the experimental number of r Bd ~ 1-53 ps. 
The Rome group [20] used a different normalisation, guided by the wish to eliminate the huge 
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uncertainty due to f Ba \ AT S /T S = (AT s /AM s ) theoi y(AM s /AM d ) theov yAM° xp T Bs . The values 
obtained by the Rome group for Ar s /r s were typically considerably lower than 0.10, which was 
partially due to different input parameters like the bottom mass. Since now AM S is known ex- 
perimentally one can abbreviate their method to Ar s /r s = (Ar s /AM s ) theory AM s exp rs s . This 
prediction assumes that no new physics effects contribute to the mass difference. This is numeri- 
cally equivalent to the use of f Ba = 230 MeV in our approach (see the passage below Eq. (l45l)). 
With that input we obtain from our analysis Ar s /T s = 0.10 ± 0.06 which is in perfect agree- 
ment with the latest update of the Rome group from this year [41]. Thus we see no discrepancy 
anymore between our predictions and those of the Rome group. 

However, our predictions have been criticised recently in [12]. The authors of [12] obtain a 
much lower central value - Ar s /T s = 0.067 ± 0.027 - and claim that this difference stems from 
their use of lattice values for the l/m& -operators, while in our approach the vacuum insertion 
approximation was used. Lattice values for the l/m& corrections can be extracted from [22] for 
the operators R , Ri and R\, but their use does not resolve the numerical discrepancy. With the 
help of one author of [12] we have traced the difference back to the omission of the radiative 
corrections contained in a\ and a 2 , when Eq. (|29l) is used to extract (R ) from lattice data on 
(Q), (Qs) and (Qs)- This is numerically equivalent to shifting B Ro from 1.1 to 1.7. If we use 
this number and f Bs = 230 MeV we obtain Ar s /T s = 0.079, which is closer to but still larger 
by 18% than the value obtained in [12]. 

Now we turn to the results in the new basis: For a direct comparison with the old operator 
basis, we first show results for the scheme characterised by m pole and z: 



^■ppolc 



Ar po1 



e,LO 



polc,s 



+ 



AI\ \ po1 



AM, 



ft 



0.0955 + 0.023.B; 



v 240 MeV / 

(0.033-B^ - 0.0065^ + 0.0055 R 

— — — \ [0.121S + 0.029-B'c,- 
x 240 MeV/ L s 

(0.0335^ - 0.0065^ + 0.005B fl 

B' Q ^Bn ( A 

12.9 + 0.5-^ + 1.7-^ ' 
B B 

hi 

0.20 + 0.02-^ + 0.44 



ps 



-1 



41.4 + 10.0- 



B 



14.4 




(48) 



• 10" 4 (49) 



Now we are in the desired situation that AY S is dominated by B and the lion's share of Ar s /AM S 
can be determined without any hadronic uncertainty! Moreover the size of the 1 /m^-corrections 
has become smaller, because the magnitude of the contribution from R is reduced by a factor 
of 3 (as anticipated from Eqs. (|25T) and (|26|) ) and the sign of this contribution has changed. 
We are left with a l/m b correction of 22% of the LO value or 28% of the NLO-value. Using 
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the new operators the a s -corrections have become smaller (22% of the LO value), too, and the 
unphysical /^-dependence has shrunk. In the case of af s = Im (rf 2 /Mf 2 ) the situation did not 
change much due to the change of the basis. Here we have no strong recommendation on what 
basis to choose. However, in the presence of new physics af s also involves Re (rf 2 /Mf 2 ) and the 
same improvements occur, as discussed in Sect. HI 

Using the non-perturbative parameters from set I we obtain the following number for Ar s : 



AI\ 



(0.081 ±0.036) ps" 



Ar, 



Ar s ■ r Bd = 0.124 ±0.056 



(50) 



The central value in the new basis is larger than the old one, while the theoretical errors have 
shrunk considerably. The numerical difference stems from uncalculated corrections of order 
a s /mb and a 2 s . As a consistency check of our change of basis one can compare the results in 
the old and the new basis neglecting all l/rrib and a s -corrections and setting B = 1 = B' s . As 
required we get in both cases the same result: Ar s /T s = 0.1497. 

For our final number we still go further. First we sum up logarithms of the form zlnz 
by switching to schemes using z defined in Eq. (l30l) . Second we calculate our results for two 
schemes of the b-quark mass, using either m& or m^ ole of Eq. (l32l) and finally average over the 
schemes. By this we obtain the main result of this paper: 



Ar, 



B, 



240 MeV , 



(0.105 ± 0.016)5 + (0.024 ± 0.004)5^ 
(0.030 ± 0.004)5^ - (0.006 ± 0.001)5% + 0.0035 K )] ps" 1 (51) 



B' B 
(9.7 ± 1.6) ±0.3-f ±0.3—^ 
B B 




10" 



+ 



Ar 5 

AM, 



B' B 
(0.08 ± 0.01) + 0.02-^ + (0.05 ± 0.01) — 
B B 

B' 

(46.2 ±4.4) + (10.6 ± 1-0)-^ 




10" 



13.2± 1.3)- 



B 



R 2 



B 



(2.5 ±0.2 



B 



B 



+ (1.2 ±0.1)^ 



10" 



Using the parameter set I, we obtain the following final numbers 

Ar., 



Ar, 



Ar s 

A Mo 



(0.096 ± 0.039) ps" 1 = 

1 s 

(2.06 ±0.57) • 10" 5 

(49.7 ±9.4) • 10" 4 

(4.2 ± 1.4) • 10" 3 = 0.24° ±0.08° 



Ar s ■ r Bd = 0.147 ±0.060 



(52) 



(53) 

(54) 
(55) 
(56) 
(57) 



The first striking feature of these numbers is the large increase for the prediction of Ar s from 
0.070 ps" 1 to 0.096 ps" 1 (about 37 %). The change of the basis is responsible for an increase of 
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about 16 %. We have shown that the previously used basis suffers from several serious drawbacks 
— most importantly in the old basis strong cancellations, which are absent in the new basis, 
occur. Next we have reduced an additional uncertainty by summing up logarithms of the form 
z In z to all orders. This theoretical improvement results in another increase of about 11%. The 
averaging over the pole and MS schemes results in an increase of about 7% compared to the 
exclusive use of the pole-scheme. Finally we also include subleading CKM- structures (as done 
in [20,21] as well) giving an increase of Ar, by about 3% compared to setting V uh to zero. In 
the case of the flavour- specific CP-asymmetry the choice of the new basis has no dramatic effect. 
If one assumes that there is no new physics in the measured value of AM, one can avoid the 
large uncertainty due to fs s by writing: 



/ Ar \ Theory 

Ar, = -r-rr -AAf Exp - = 0.088 ±0.017ps~ 



(58) 



^ = AT s -r Bd = 0.127 ±0.024. (59) 

J- s 



This smaller value is numerically equivalent to using f Bs = 230 MeV in Eq. (l5D) . 
For completeness we also present the numbers with the parameter set II: 



Ar, 



AM, 



(0.106 ±0.032) ps" 1 
(2.06 ±0.57) • 10" 5 
(51.9 ±9.8) ■ 10~ 4 



Ar, 



AVs-tb, = 0.162 ±0.049 



(60) 
(61) 
(62) 



The above errors in Ar, and AM, have to be taken with some care, since we were not using our 
conservative error estimate but the preliminary values from [24] . 

In the following table the individual sources of uncertainties in Ar, — using the parameter 
set I — are listed in detail: 



3 Numerical predictions 



17 



Input 










old, pole, z 


new, pole, z 


new, average, z 


fs a 


i +0.361 
1 -0.306 


-, +0.361 

-■--0.306 


i +0.361 

-■-—0.306 


Bi 


1 ±0.002 


-i +0.070 

— U.U / 1 


1 ±0.066 


±J0 '3. 
— z,o 


1 ± 0.167 


1 ± 0.035 


1 ± 0.031 


it-?. 


1± 0.235 


1± 0.203 


1± 0.157 


Br 


1 ±0.140 


1± 0.036 


1 ±0.030 


/ii with m b /2 < Hi < 2m b 


i +U.248 
- L -0.521 


1 +U.111 

x -0.272 


i +0.UY4 
-■--0.200 


v cb 


_ in / i/iu 
-i -r-u.u^ty 

-■--0.048 


, III /HQ 

-■--0.048 


1 ± 0.049 


z 


-. 4-11 044 
i -0.046 


-r 4-0 040 

-■--0.042 


1 ± 0.019 


m b 


i +0.043 
i -0.042 


-i +0.036 
-■--0.035 


-1+0.010 
-■--0.009 


a s 


1 +0.014 

-■--0.013 


1 ±0.003 


1 ±0.001 


m s 


1 ±0.010 


1 ±0.012 


1 ±0.010 


7 


-I +U.UU5 
-■--0.016 


-i +0.005 
-■--0.015 


i +0.005 
-■--0.014 


\V ub /V cb \ 


1 ±0.006 


1 ±0.006 


1 ±0.005 




1± 0.607 


1 ± 0.450 


1 ± 0.405 


now 
m b 


1-0.368 


1-0.158 


1-0.112 


RS 


1 ±0.133 


1 ±0.065 


1 ±0.066 



(63) 



The same result is visualised in figure [3] In the case of Ar s the by far largest uncertainty 
stems from the error on fs s ■ Here a considerable improvement from the non-perturbative side 
is mandatory. The dependence on the decay constant is of course not affected by the change of 
the operator basis. The second most important uncertainty comes from the l/m^-operator R 2 . 
This operator has up to now only been estimated in the naive vacuum insertion approximation. 
Any non-perturbative investigation would be very helpful. Number three in the error hit list is 
the unphysical /^-dependence. Using the old operator basis the corresponding error was huge, it 
was drastically reduced by changing to the new basis and by including also the MS-scheme for 
the b-quark mass. Any further improvement requires a cumbersome NNLO calculation, which 
might be worthwhile if progress on the non-perturbative side for /s a and R 2 is achieved. Number 
four is again a non-perturbative parameter - now the bag parameter of the operator Q. In the old 
operator basis the corresponding uncertainty stemmed from Bs and was larger by a factor of 2.5. 
The dependence on V cb results in a relative error of about 5% for both the old basis and the new 
basis. All remaining uncertainties are at most 3%. 

Using our conservative estimates and adding all errors quadratically (after symmetrising them) 
we arrive at a reduction of the overall theoretical error due to the introduction of the new basis 
from ±61% to ±41%, where the last number is completely dominated by the decay constant. 
If one neglects the dependence on /s a the overall theoretical error goes down from ±51% to 
±23%. 

In the table in Eq. (1631) we also show the dependence on the b-quark mass we are using in the 
l/m fe -corrections, m^ ow . This dependence can be viewed as a measure of the overall size of the 
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Figure 3: Uncertainty budget for the theory prediction of AT S . The largest uncertainties stem 
from f Ba , the renormalisation scale /ii of the AB = 1 operators and the bag parameter of 
the l/mt-suppressed operator R 2 . The transparent segment of the right pie chart shows the 
improvement with respect to the old result on the left. 



l/m & -corrections. The use of the new basis results in a strong reduction of the corresponding 
uncertainty, from 37% to 11%. And finally we compare the two renormalisation schemes (RS) 
we are using for the b-quark mass. Here we have again muss less uncertainty in the new operator 
basis. To avoid a double counting of the errors we did not include the last two rows of table (1631) 
in the total error. 

Investigating the case of AT S /AM S the improvement due to our new basis is more substantial, 
since here the dependence on f Bs cancels: 
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Input 


AT JAM. 

Sj — S 


AT J AM, 

Sj — S 


a? 

is 




old, pole, z 


new, average, z 


new, average, z 


1 - > \ 


i +U.U74 
1 -0.064 


1 ± 00^ 

X _l_ U.UUu 


i +0.006 
-■--0.005 


^2,3 


i + o 1 67 


1 ± 031 


1 + 004 

X _L W.WWt: 


R - 


1 _l n ooc 
1 ZL U.ZiOU 


i ±n lev 
i nz u. i o i 


i _l n r>\ 

1 ZL U.UZO^il/jJ 


R„ 


i _l n i a n 

1 ZL U . X'rfcU 


i _l n non 
1 ZL u.uou 


X ZL U.UXXIXL3I 


1 1 -, xx/ 1 1 n m 7 / <C 1 1 -, <C m , 

Ll\ Willi III}} ^ Ub\ ^ ^lll\y 


1 +U.194 
1 -0.495 


-i +0.027 
-■--0.154 


1+0.152 

-■--o.ioi 


V V, 

'CO 


1 +0 000 

X _l — \J t\J\J\J 


1+0 000 

X _l — w.www 


1 + n 000 

X _l — \J •\J\J\J 


z 


-I +0.044 
-■--0.046 


1 ±0.019 


i +0.094 
-■--0.092 


m b 


i +U.043 
-■--0.042 


i +0.01U 

-■--0.009 


i +0.03Y 
-■--0.036 


m t 


1± 0.018 


1± 0.018 


1 + 0.018 


a s 


1± 0.012 


1± 0.001 


1± 0.007 


m s 


1± 0.010 


1± 0.010 


1 ±0.001 


7 


i +U.0U1 
-■--0.003 


i +0.00U 

-■--o.ooi 


1+0.144 
-■--0.081 


\v ub /v cb \ 


1 ±0.001 


1± 0.001 


1+0.194 
-■--0.196 




1 ± 0.480 


1± 0.189 


1 ± 0.279 


now 

m v b 


1-0.368 


1-0.112 


+0.016 


RS 


1 ±0.136 


1 ±0.069 


1± 0.004 



In the case of Ar s /AM S the use of the new operator basis leads to a reduction of the total 
error from 48% to 19%! The dominant error is now due to the bag parameter B R2 , followed by 
the /ii -dependence. The remaining uncertainties are at most 3%. In the case of af s the situation 
is quite different. Here the dominant uncertainty stems from V ub , followed by the dependences 
on fii, 7 and z. Moreover the l/m^-corrections play a minor role here — as can be read off from 
the error due to the variation of m^ ow . 



3.4 AM d , AT d and af s within the SM 



Here we give updated numbers for the mixing parameters of the B d system. The CKM elements 
governing B d — B d mixing appear in the combinations Xf = V* d V ib for % = u,c, t. The bag 
parameters multiplying f Bd below refer to B d mesons and are different from those in the B s sys- 
tem. However, no non-perturbative computation has shown any numerically relevant deviation 
of B Bd /B Bs from 1. 

Updating AM d to m t (m t ) = 163.8 ± 2.0 GeV gives 



AM d = (0.53 ± 0.02) ps 



-l 



td\ 



B 



0.0082 



200 MeV/ 0.85 



While in the B s system the values of 7 and \ V ub \ in Eq. (|37l) play a minor role, their uncertainties 
are an issue for AT d and at. The master formulae are [21] 



AT d 
AM, 



-10" 



c + aRe _| + 6Re _|_ 



(64) 
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Figure 4: Uncertainty budget for Ar s /AM S . See Fig.[3]for explanations. The ratio Ar s /AM S 
does not depend on f Bs and the progress due to the new operator basis is more substantial than 
in Ar s . 
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The coefficients 
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(66) 



are independent of CKM elements because of Mf 2 oc Xf 2 . In our new operator basis these 
coefficients read 



,+0.09\ -B5 



0.3Hao; + 0.271^ 



0.15\ ^i? 



5 



0.08 ± 0.03 + (0.02 ± 0.01) ^ + (0.04^ J?) % 

B \ ' B 



-46.1 ± 6.6 - (10.5 ± 1.3) ^ + (8.7^o) B 
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With the hadronic parameters of Set I in Eq. (l39l) one finds 

b = 0.2 ±0.1, 



10.51};?, 



,o 0+12.7 



(67) 



It is convenient to express X^/Xf in Eqs. (|64|) and (1651) in terms of the angle (3 
of the unitarity triangle and the length R t = \Xf/X^\ of the adjacent side [21]: 



arg(-Af/A 



Re^ 
X? 

\ d 
A? 



cos (3 
sin /3 



i. 



Re 



A 



J 2 



Im 



Af 



cos(2/3) 2 cos/3 + l 



^ 2 



sin(2/3) 



+ 2 



sin /3 



(68) 



Clearly the terms involving Xf^/Xf 2 in Eqs. (|64l) and (1651) are numerically irrelevant in view of 
the smallness of b. Moreover, in the preferred region of the Standard Model fit of the unitarity 
triangle one has cos (3 ~ R t , so that Re Af / Xf is suppressed. Setting a and b to zero in Eq. (l64l) 
reproduces Ar d /AM d within 2% [21] and AY d / AM d is essentially free of CKM uncertainties. 
Inserting Eqs. (|67]) and ([68]) into Eqs. ([64]) and ([651) yields 



AM, 



53.3l£ 4 + (10.3+1 



8 

-1.7 



1 



cos(/3) 



+ (0.2 ±0.1] 



'cos(/3) cos(2/3)' 



10" 



(69) 



4 



sin /3 



+ (0.2 ±0.1 



sin(2/3)' 

— w 



(70) 



Next we insert the numerical values for (3 and R t from [2]. Since we are interested in testing the 
hypothesis of new physics in B s — B s mixing, we take values for (3 and R t obtained prior to the 
measurement of AM S . With (3 = 23° ± 2° and R t = 0.86 ±0.11, which correspond to a CL of 
2a, one finds 



AM, 



52.612:5) • io- 4 , 



-4.8 



+1.0 
-1.2 



10" 



(71) 



Thus these predictions allow for new physics in AM S , but assume that all other quantities 
entering the standard fit of the unitarity triangle in [2] are as in the Standard Model. Using 
0.507 ± 0.004 ps" 1 and r^ xp = 1.530 ± 0.009 we find from Eq. (HQ): 



AMf p 



Ar n 



Ar 

AMT P 
AM, d 



(26.7 



+5.8' 
-6.5 



10" 4 ps" 



ar 



(40.9lg;g) • 10" 4 . (72) 



The result in Eq. (1721) is consistent with our prediction in [21], but the central value is substan- 
tially higher. This is not solely caused by our new operator basis, but also by the use of a different 
renormalisation scheme. In both [21] and this work we average over two schemes, but in one 
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of the schemes used in [21] the zln z terms are not summed to all orders. Note that the quoted 
error of af s in [21] corresponds to the la ranges of (3 and R t , while in Eq. dTTI ) more conser- 
vative 2a intervals have been used. The ranges in Eq. (1711) imply for the CP-violating phase 

d = arg(-M 1 d 2 /rf 2 ): 

<p d = -0.091™ = -5.2°^° • (73) 



4 Constraining new physics with B s — B s mixing 

In this section we investigate effects of new physics contributions to the £> s -mixing parameters. 
New physics can change the magnitude and the phase of M* 2 . We parameterise its effect (simi- 
larly to [2,42]) by 

M-{ 2 = Mf 2 M ' s • A s , A s = |A s |e^. (74) 

The relationship to the parameters used in [2, 14] is 

A s = r 2 s e 2id *. 

We find it more transparent to plot Im A s vs. Re A s than to plot 26 s vs. r 2 . Our plots are similar 
to Fig. 1 of [14], which displays sin(20 s ) vs. cos(2# s ), but also include the information on | A s | = 
r s . Finally T s 12 stems from CKM-favoured tree decays and one can safely set T\ 2 = T™' s . 

4.1 Ar s , £T S /AM S and a| beyond the SM 

One easily finds: 

AM S = AM S SM |A S | = (19.30 ± 6.74) ps _1 • |A S | (75) 
Ar s = 2|r^ 2 | cos (<^ M + A ) = (0.096 ± 0.039) ps" 1 ■ cos (<^ M + A ) (76) 

Ar s |r? 2 | cos(0sm + 0^ 



= v° iv = (4.97 ± 0.94) -10- 3 v ° , ° y (77) 



cos ( 


0f* + 




|A S | 



|Pf 2 | sin(0f M + 



3 Sinl 











< = 7 ± ^£h a i = (4.97 ± 0.94) -10- d , . , J (78) 

|M 12 ' | |A S | v 7 1 v 

with (cf. Eq. d57j) <^ M = (4.2 ± 1.4) • 10~ 3 (79) 

Here the numerical values correspond to our results from parameter set I in Eqs. (I54l457l) . In 
the case of af s there is a major difference to the SM case of Sect. 13 .3 1 which only involves 
Im(rj 2 /Mf 2 ): in the presence of new physics af s is dominated by Re(r^ 2 /M{ i 2 ) as long as 
\(f)f | > 0f M . Thus the prediction in Eq. (1781) profits from the improvements due to our new 
operator basis — just as the prediction of Ar s in Eq. (1771) . From Eq. (1781 one also verifies the 
enormous sensitivity of a| to new physics, since it exceeds its SM value by a factor of 250 for 
4>f = tc/2. We have plotted af s vs. A for the old and the new bases in Fig. [51 
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Figure 5: af s as a function of the new phase </>f from Eq. (|78T) for the range — n < (fif < ir. 
The thick blue lines show the prediction in the new basis, while thin red lines correspond to the 
old operator basis. The solid lines display the central values of our predictions and the dashed 
lines show the uncertainties, which are much larger for the old result. The standard model value 
af s (0f = 0) = 2.1 • 10" 5 is too close to zero to be visible in the plot. 



4.2 Basic observables 

In this section we summarise the observables which constrain |A g | and <pf. These constraints 
are illustrated in Fig.[6]for hypothetical measurements. 

1. The mass difference AM, determines |A S | through Eq. (|75T) . The accuracy of |A S | ex- 
tracted from AM, is limited by the precision of a lattice computation. This is not the case for the 
other quantities discussed in this section. 

Alternatively one can confront the experimental ratio AM d / AM S with theory. This has the 
advantage that the ratio of the hadronic matrix elements involved can be predicted with a smaller 
error, of order 5%. However, then the parameter of R t of the unitarity triangle entering AM d 
must be taken from measurements which are insensitive to new physics (or at least insensitive 
to new physics in B s — B s mixing), e.g. through determinations of the CKM angle 7 from tree- 
level B decays (cf. the discussion after Eq. (1701). At present this method leads to comparable 
uncertainties in the extracted |A S | as the direct determination from AM S . (Further flavour-blind 
new physics cancels from AMd/ AM S .) In the following analyses we do not use AMd/ AM S . 

2. The lifetime measurement in an untagged b — > ccs decay B s — > fcp, where fcp is a CP 
eigenstate, determines Ar s cos(0f - 2f3 s ) = \ AT S cos(0f - 2f3 s ) \ [43,44]. Consider a CP-even 
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final state fcp+ like DfD s . The time-dependent decay rate reads 

r[ ( ^ -+ W] oc 1 + cos( f -W'h-™ + i-«»(^-2A) e _ rjrt 

cosh — cos(0 s — 2p s ) sinn 



-r.t 



(80) 



2 ' s/ 2 

and the (time-independent) overall normalisation is related to the branching fraction [44]. Here 



Ps = -arg = 0.020 ±0.005 = 1.1° ±0.3°. (81) 




That is, — P s is the analogue of the angle (3 of the unitarity triangle, which governs the mixing- 
induced CP asymmetry in Bd — > J/iftKs, in the B s system. For (3 S different sign conventions 
are used in the literature, we chose the one of [6] which satisfies (3 S > 0. 

For example within the Standard Model (and neglecting the tiny [3 S ) the lifetime measured in 
( £>s — ► D+D~ equals T S L = T s + AT/2, because only the short-lived CP-even mass eigenstate B L 
can decay into D+D~. By using the theory relation l/re d = Yd = (1.00±0.01)r s one then finds 
Ar s . For (f)f 7^ 0, however, the mass eigenstates are no more CP eigenstates and both of them 

can decay to a CP eigenstate, as can be easily verified from Eq. (|80l) . From r[ B s — ► fcp+, t] one 
can extract | Y s | , | AT S | , | cos(0f ) | and the overall normalisation, if the statistics is high enough to 
separate the two exponentials. If the measured Y[ { B } S — > fcp+, t] is fitted to a single exponential 
exp[— Tf t], the measured rate is [44,45] 



1 + cos(<# - 2ft))/r L + (1 - cos(0f - 2f3 s ))/T H 

(Ar s ) 2 ^ 



Tf (i + cos(0f _ 2 /3 s ))/rl + (i - cos(0f - 2/3.))/r* B (82) 



= Y s + Ar s cos(0f - 2(3 S ) + o 

= T s + 2 \T S 12 \ cos(0f + S SM ) cos(0f - 2/3.) + • (83) 

For a CP-odd final state one has to interchange Y l and Yh in Eqs. (f8~0l) and (f82l) and to flip the sign 
of cos(0f — 2/3 s ) in Eqs. (|80l) and (1831 . From Eq. (1831) it is clear that the lifetime measurement 
determines [43,44] 

Ar s cos(0f) = 2|ry C os 2 (0f), 

if the small phases 0f M and (3 S are neglected. Thus one can find | cos^f |, which determines cf>f 
with a four-fold ambiguityfl We stress that (since sign Ar s = sign cos(0^)) the lifetime method 
gives no information on the sign of Ar s and experimental results should be quoted for |Ar a | 
rather than Ar\. 



*If one keeps <j>^ M and j3 s non-zero, one solution for cf>f is related to the other three by <j)^ — > <pf + tt, <p^ 

2f3 s - 0™ _ A and(/) A ^ Ws _ SM _ A +7r _ 
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Eq. (f8~2l) assumes that detection efficiencies are constant over the decay time. Since this is not 
the case in real experiments, we strongly recommend to perform a three-parameter fit to 2|r 12 |, 
| cos(0f )| and the overall normalisation (with |r s | fixed to (1.00 ± 0.01)) to Eq. (1801) . 

With the advent of the precise measurement of AM S [10] one will rather exploit | Ar s |/AM S 
to constrain A s than |Ar s | itself, which suffers from much larger hadronic uncertainties. From 
Eq. (1771) one infers that | Ar s |/AM S defines two circles in the complex A s plane which touch the 
y-axis at the origin. 

3. The angular analysis of an untagged b — ► ccs decay { B } S — ► VV, where VV is a superposi- 
tion of CP eigenstates with vector mesons V, V, not only determines Ar s cos(</>f — 2f3 s ), but also 
contains information on sm((pf— 2(3 S ) through a CP-odd interference term. Here the golden mode 
is certainly B s — ► J/ip(j), but also final states with higher ip resonances and B s — > D* + D*~ can 
be studied. The determination of 4>f from the CP-odd interference term in untagged samples 
involves a four-fold ambiguity. It could be reduced to a two-fold ambiguity if the signs of cos 5\ 
and cos 5 2 were determined, where 8\ and 5 2 are the strong phases involved [44, 46] . These two 
solution are related by (j)f <-> ± it. If one relaxes the assumptions on cos S 1 and cos 5 2 , one is 
back to the same four-fold ambiguity as in item 2. 

4. The branching fraction Br( ( B } s — ► D^+D^,*^) approximates the width difference AFcp 
between the two CP eigenstates of the B s system [44]. Irrespective of any new physics in M( 2 
one always has Ar CP = 2|r* 2 |, so no constraint on our new physics parameter A s is gained. 
Yet the ratio of Ar s cos(0f — 2(5 S ) and AFcp could cleanly determine cos(0f ) cos(0f — 2(5 S ). 
However, Br(B s — > £)(*)+£)(*)-) only equals Ar CP in the poorly tested simultaneous limit of 
an infinitely heavy charm quark with small-velocity [47] and an infinite number of colours [48]. 
In order to test this limit one needs to measure the CP-odd and CP-even fractions of all b — > ccs 
decays [44]. Until this has been done nothing can be inferred from Bt(b ] s — > D^ + D^~), 
in particular this quantity neither gives an upper bound (since other CP-even b — > ccs modes 
can be relevant) nor a lower bound (since other CP-odd b — > ccs modes can be relevant and the 
£)(*)+£)(*)- fl na i state has a CP-odd component) on Ar CP . We strongly discourage from the 
inclusion of Br( B s -> D^ + D^~) in averages with Ar s determined from clean methods. 

5. af s can be measured from untagged flavour- specific B s decays, typically from the number 
of positively and negatively charged leptons in semileptonic decays. Observing further the time 
evolution of these untagged ( B } S — > ( v\ decays (see e.g. [8]), 



will have two advantages: one can use the oscillatory behaviour to control fake effects from 
experimental detection asymmetries (which are constant in time) and to separate the B s and B d 
samples through AM S ^ AM d . The constraint from af s on A s is given in Eq. (|78l) . It defines a 
circle in the complex A s plane which touches the x-axis at the origin. The constraint from a| s 
on A s only has a two-fold ambiguity (related to 4> s «-> ix — 4> s ) and discriminates between the 
solutions in the upper and lower half-plane in Fig. [6l 

6. The time dependence of the tagged decay B s — > J/ip(j) permits the determination of the 
mixing-induced CP asymmetries A™p(B s — > (J/4>(fi)cp±)- The angular analysis separates the 



rfgj -> x-e+u e ,t] - r[ ( £?j -> x+e-v t ,t] 
TfEl -> x-e+v t ,t] + rfBl -> x+e-v t ,t] 




cosh (AT s t/2) 



cos (AM, t) 



(84) 
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CP-odd P-wave component from the CP-even S-wave and D-wave. The time-dependent CP 
asymmetry is (in the notation of [6,44]): 

r(B?(t) -> f) - T(B° a (t) -> f) A™sm(AM s t) 

T(B°(t) - /) + r(5°(t) - /) cosh(Ar,t/2) + A Ar sinh(AlV/2) ' 

One finds 4>f — 2(3 S through 

A^(B S -> (J/^)cp±) = ±sin(«^-2ft), A Ar = T cos(0f -2&) (86) 

with the same two-fold ambiguity as from af s in item 5. Combining Eqs. (|75|) and (1781) with 
Eq. (l86l) and neglecting the tiny contributions of and (3 S one verifies the correlation between 
af s and A™p x (£> s — >■ (J/ip4>)cp±) derived in [12, 13]. In fact such correlations can be found 
between any three of the observables discussed above, because the B s —B s mixing only involves 
the two parameters |A S | and <p s . 

An important remark here concerns the decay B s — > K + K~, as one might be tempted to use 
the lifetime measured in B s — > K + K~ to determine r s + |Ar s /2|. While K + K~ is CP even, the 
decay is penguin-dominated and as such sensitive to the same kind of new physics which may 
be responsible for the experimental anomaly seen in penguin-dominated Bd decays [3]. Thus 
information from B s — > K + K~ should under no circumstances be included in any averages 
with the measurements discussed above. Instead one should confront the lifetime measured in 
this mode with the one obtained from B s — ► (J/ijj(f))cp+ to probe new physics in b — > s penguin 
decays. 

For a visualisation of the bounds from Eqs. (17514781 ) in the complex A s -plane we consider 
now the hypothetical case of |A S | = 0.9 and 4>f = — n /A. Suppose one would measure these 
central values: 

AM S = 17 A ps~\ Ar s = 0.068 ps^ 1 , (87) 

— ^- = 3.91 • 10~ 3 , at = -3.89- 10~ 3 . (88) 

AM S ts 

Moreover we assume the following theoretical and experimental uncertainties: AM S : ±15%, 
A(f) s : ±20%, AVs/AMs : ±15%, af s : ±20%. The regions in the A s -plane bounded for these 
hypothetical measurements are shown in figure [6[ 

The constraints from CP-conserving quantities are symmetric to the Im(A s )-axis, The bound 
from AM, simply gives a circle with the origin (0,0) and the radius |A S |. In the measurement 
of Ar s we have assumed that the data are fitted to the correct formula Eq. (l80l) and |Ar s | and 
| cos (0 S — 2{3 S ) \ have been determined as discussed above in item 2. In practice the extracted 
| Ar s | and | cos((fif — 2(3 S ) \ are strongly correlated and mainly | Ar s | | cos(0f — 2(3 S ) \ is determined 
(see Eq. (|8~3l and [44]). The constraint from the hadronically cleaner ratio |Ar s |/AM s are 
two circles which touch the y-axis in the origin. If one fully includes the correlation between 
|Ar s | and | cos(0f — 2(3 S ) \ one will rather find constraints which roughly correspond to a fixed 
|Ar s cos(0f — 2[3 S )\/ AM S . The corresponding curves are a bit more eccentric than the circles 
from |Ar s |/AM s . 
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Figure 6: Illustration of the bounds in the complex A s -plane for | A a | = 0.9 and <pf = — tt/4. We 
assume the following overall uncertainties: AM S (red or dark-grey) : ±15%, Ar s /AM S (yellow 
or light-grey): ±15%, af s (light-blue or grey) : ±20% and <j>f (solid lines) : ±20%. 



If one plots the bounds from |Ar s | (or |Ar s cos(0^—2/3 s )|) alone, one finds four rays starting 
from the origin. The experimental information in this is redundant, as it is fully contained in 
the constraints from AM S and |Ar s |/AM s . For the theory uncertainties, however, this is not 
true: if (as current data do) AM S prefers a small value of fs s , while Ar s prefers a large fs s , 
the combined constraint from AM S and |Ar s | will exclude a region of the A s plane which is 
allowed by the ratio | Ar s |/AM S , from which f Bs drops out. 

The measurement of af s yields a circle touching the x-axis in the origin, in particular it 
reduces the four-fold ambiguity in the extracted value of A s to a two-fold one. The extraction of 
4>f — 2f3 s from the angular analysis in B s — > J/ip<f) (as discussed in item 3) also yields four rays 
starting from the origin (corresponding to the same value of | cos(0f — 2(5 S ) |), if no assumptions 
on the signs of cos 5\ and cos 82 are made. Finally, the measurement of Aq^(B s — > (J/ip4>)cp±) 
will select two out of these four rays, discriminating between <fif — 2f3 s > and <fif — 2f3 s < 0. 

4.3 Current experimental constraints on A s 

In this section we turn to the real world and discuss the current experimental constraints on the 
complex A s -plane. In view of the experimental errors we set (3 S to zero and identify <p s with (fif. 

The mass difference AM, is now known very precisely [10], see Eq. ©. For the remaining 
mixing parameters in the £> s -system only weak experimental constraints are available. The only 
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available experimental analysis of | Ar s | with the correct implementation of the phase <p s is from 
the D0 collaboration, their analysis in [49] was recently updated in [50] using lfb -1 of data. 
Setting the value of the mixing phase S to zero (Standard Model scenario) they obtain [50] 

Ar s = O.^iO.OS^lSS^ps" 1 . (89) 

Allowing for a non-zero value of the mixing phase S they get 

Ar s = 0.17±0.09 (stat) ± 0.03 (syst) ps- 1 

and cj) s = -0.79 ± 0.56 (stat) ± 0.01 (syst) (90) 

or Ar s = -O.lTiO.Og^iO.OS^ps- 1 

and (f) s = -0.79 ± 0.56 (stat) ± 0.01 (syst) + vr . (91) 



As expected from Eq. (183]) the values for | Ar s cos <p s \ found from Eqs. (1891 ) and (1911 ) are roughly 
equal to Ar s in Eq. (l89l . The quoted results in Eqs. (l90l) and (19)1 assume that the signs of cos Si 
and cos 5 2 agree with the results found with naive factorisation. With this assumption the other 
two solutions for S (which have opposite signs to those in Eqs. (|90l) and (19TI) ) are excluded. 
Strategies to check this theoretical input are discussed in [44]. 

The semileptonic CP asymmetry = a| s in the B s system has been determined directly 
in [51] and was found to be 

^direct = (24.5 ± 19.3 (stat) ± 3.5 (syst) ) • 10~ 3 . (92) 

Moreover the semileptonic CP asymmetry can be extracted from the same sign dimuon asym- 
metry that was measured in [52] as 

a al = (-5.3 ± 2.5 (stat) ± 1.8 (syst) ) ■ lO" 3 (93) 

in a data sample containing both Bd and B s mesons. While the composition of the sample is 
known, no determination of the initial state on an event-by-event basis was possible. Updating 
the numbers in [14,53] one sees that the measurement in Eq. (|93T ) determines the combination 

a sl = (0.582 ± 0.030) a d sl + (0.418 ± 0.047) a s sl . (94) 

In [14, 53] the experimental bound for from B factories was used to extract a bound on 



from Eq.(l93l and Eq.(l94l). The huge experimental uncertainty in then inflicts a large error on 
the value of inferred from Eqs. (|93T ) and (|94| ). 

Here we pursue a different strategy and use the much more precise theoretical Standard 
Model value for in Eq. dTD) . In the search for new physics this is permissible: if the re- 
sulting constraint on A s departs from the Standard Model value A s = 1, this will then imply 
new physics in either afj or a^. Moreover, the current precision in the unitarity triangle already 
substantially limits the room for new physics in [2]. 

Using 4 = - (0.48^;{°) • 10~ 3 of Eq. (ED and further Eqs. ([93]) and ([94]) we obtain the nice 
bound 

a s,dimuon = ( _ 12 Q ± g Q(stat) ± 4 5(syst)) . 10 _ 3 (95) 
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Figure 7: Current experimental bounds in the complex A s -plane. The bound from AM S is given 
by the red (dark-grey) annulus around the origin. The bound from |Ar s |/AM s is given by the 
yellow (light-grey) region and the bound from af s is given by the light-blue (grey) region. The 
angle <pf can be extracted from | Ar s | (solid lines) with a four-fold ambiguity — each of the four 
regions is bounded by a solid ray and the x-axis — or from the angular analysis in B s — > J/\P0 
(dashed line). This constraint also has a four-fold ambiguity if no assumptions on the strong 
phases 5\ and <5 2 are made. The dashed lines limit the region corresponding to the solution in 
Eq. (l90l ). The Standard Model case corresponds to A s = 1. The current experimental situation 
shows a small deviation, which may become significant, if the experimental uncertainties in Ar s , 
a s sX and 4> s will go down in near future. 



Combining this number with the one from the direct determination [51] in Eq. (1921) we get our 
final experimental number for the semileptonic CP asymmetry: 

a s sl = (-8.8 ± 5.7 (stat) ± 4.5 (syst) ) • 1(T 3 . (96) 

Adding statistical and systematic error in quadrature gives 

a' A = (-8.8 ± 7.3) ■ 1(T 3 . (97) 

In Fig. CD) we display all bounds in the complex A s -plane including all experimental and theo- 
retical uncertainties. 

The combined analysis of AM S , <f> s , |Ar s |/AM s and in Fig. [7] shows some hints for 
deviations from the Standard Model. To analyse them further we ignore discrete ambiguities 
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and focus on the solution in the fourth quadrant which is closest to the Standard Model solution 
A s = 1. We further do not perform a complete statistical analysis with proper inclusion of all 
correlations and for simplicity add statistical and systematic errors in quadrature. First we note 
from Eq. (1751) that Eq. © implies 



| A, | = 0.92 ± 0.32 (th) ± 0.01 (exp) (98) 



while Eqs. C3) and ([90]) lead to 



COS( 



= 1.93±0.37 (th) ±l.l (exp) . (99) 

Eqs. (|98"1) and (|99~1) are consistent with A s = 1, but prefer |A a | < 1. 

Second we observe that both the angular distribution in B s — > J/ip(J) giving Eq. (l9~TT) and 
in Eq. (1971) point towards a non-zero <p s . Both analyses involve sin S , the two values inferred are 

sin0 s = — 0.71^0 27 from the angular analysis, Eq. (|9~T1), (100) 

sin 6. . . 

— r-y = -1.77±0.33 ( th) ± 1.47 (eX p) from a sl in Eq. ([97]). (101) 

I s I 

In Eq. (11011) we have profited from our improved theory prediction in Eq. ([78]) . For |A S | = 1 the 
two numbers combine to 

sin0 s = -0.77 ± 0.02 (th) ± 0.36 ( exp). (102) 
Relaxing | A s | to its minimal value allowed by Eq. (|98"1) . | A s | = 0.59, changes this result to 

sin0 s = -0.76 ± 0.03 (th) ± 0.34 (exp) . (103) 

Either Eq. (11021) or Eq. (11031) alone imply a deviation from (fi s = by 2.1a, but Ar s in Eq. (l9~TT) 
pulls in the opposite direction, preferring large values of | cos <p s \ through Eq. d76l) . Despite of its 
large error Ar s already gives a powerful lower bound | cos0 s | > 0.55 (so that | sin0 s | < 0.84) 
at the la level because of its large central value in Eq. (|9~TI) . This can be clearly seen from Fig. IT] 
However, Ar s is consistent with cos0 s = at the 1.8cr level and clearly has no impact on the 
small <p s region, which is the relevant region to assess the significance of Eq. (11021) in the search 
for new physics. 

In conclusion we find that the data are best fit for <p s around —0.88 corresponding to sin <p s = 
—0.77, if | A s | = 1. The constraint from |AT a | is less compelling, but slightly prefers | A s | < 
and disfavours too large values of | sin0 s |. The discrepancy between data and the Standard 
Model is around 2a, which is not statistically significant yet. If our results are used to constrain 
models of new physics one should bear in mind that we have only discussed the solution in the 
fourth quadrant of the complex A s plane here. 
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5 A road map for B s — B s mixing 

Clearly the best way to establish new physics from B s — B s mixing is a combination of all ob- 
servables following the line of Sect. 14.21 In particular it has to be stressed that Aq P x (B s — > 
(J/i[)(j))cp±) and af s are not substitutes for each other, but rather give complementary informa- 
tion on the complex A s plane because of their different dependence on |Mf 2 |. With the new 
operator basis presented in this paper it will be possible to determine A s solely from measure- 
ments which involve hadronic quantities only in numerically sub-dominant terms. To this end 
any experimental progress on |Ar a |, a|, the angular distributions of both untagged and tagged 
B s — > J '/ip(j> decays (with the tagged analysis giving access to A™p(B s — > (J/ip(f>)cp±)) an d 
possibly of other b — > ccs decays of the B s meson is highly desirable. Regardless of whether 
sin cf) s turns out to be zero or not it is important to measure the sign of Ar s . Methods for this 
are discussed in [44]. Probably the most promising way to determine sign Ar s = sign cos(0 s ) 
is the study of B s — ► J/ipK + K~ with a scan of the invariant mass of the (K + , K~) pair around 
the cf) peak to determine sign cos Si >2 . 

Clearly the analysis of the precise measurement of AM, needs a better determination of 
f%B, Since any new physics discovery from a quantity involving lattice QCD will be met 
with scepticism by the scientific community, the lattice collaborations might want to consider to 
switch to blind analyses in the future. The predictions of both Ar s / AM S and af s involve the ratio 
B's/B in a numerically sub-dominant term. It may be worthwhile to address this ratio directly 
in lattice computations, because some systematic effects could drop out from the ratio of the two 
matrix elements. 

The quantities discussed in this paper will also profit from higher-order calculations of the 
short-distance QCD parts. In particular corrections of order a s /m b should be computed to permit 
a meaningful use of l/m& bag factors computed with lattice QCD or QCD sum rules. A further 
reduction of the dependence on the renormalisation scale fii requires the cumbersome calculation 
of 0(a 2 s ) corrections. Finally, the reduction of the 1/mb corrections with the help of our new 
operator basis can only be fully appreciated, if the size of the 1/ml terms is indeed small. We 
have estimated these corrections and indeed found no unnatural enhancement over their natural 
size. 

6 Summary 

In this letter we have improved the theoretical accuracy of the mixing quantity rj 2 , q = d,s, 
by summing the logarithmic terms a™zln n z, z = ml/ml to all orders n = 1,2, .. . and by 
introducing a new operator basis, which trades the traditionally used operator Qs of Eq. © for 
Qs defined in Eq. (TT6T ). In the new operator basis the coefficient of the 1 /m^-operator R is 
colour-suppressed. We have found that all previously noted pathologies in the sizes of the l/m b 
and a s corrections were artifacts of the old operator basis. Still, one could achieve the same 
accuracy with the use of the old basis, if one i) used the coefficients with resummed In z terms, 
ii) added the term of order N c a s /mb which drops from the NLO results of [19-21] when Qs is 
eliminated for R and iii) fully takes the numerical correlation between B and B s into account. 



32 



Theoretical update of B s — B s mixing 



This numerical correlation stems from the smallness of the matrix element (Qs) ■ It is most easily 
implemented by expressing either B or B s in terms of B s , which is essentially equivalent to our 
approach. 

Our improvements are most relevant for ReT^/M^, which enters both AT q /AM q and new 
physics scenarios of af s . In particular, hadronic quantities now appear in these quantities in 
numerically sub-dominant terms only. We have then discussed how experimental information on 
|Ar s |, af s , (f) s from the angular distribution of ( B S — > J/tp(J) an d A™p(B s — ► (J /V0)cp±) can 
be efficiently combined to constrain the complex parameter A s , which quantifies new physics in 
B s — B s mixing. 

Armed with our more precise formulae we have analysed the combined impact of the D0 
analyses of the dimuon asymmetry and of the angular distribution in the decay { b\ — > J / ip<j). 
Here we have assumed that cpd is free of new physics contributions. This is plausible in view 
of the constraints on (fid from global fits to the unitarity triangle [2]. Scanning conservatively 
over theory uncertainties, we find that <j) s deviates from its Standard Model value by 2 standard 
deviations. 
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